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Abstract 

We present a new cascade-type microscopic simulation of Nucleus-Nucleus Colli- 
sions at RHIC energies. The basic elements are partons (quarks and gluons) moving in 
8iV-dimensional phase space according to Poincare-covariant dynamics. The parton- 
parton scattering cross sections used in the model are computed within perturbative 
QCD in the tree-level approximation. The Q 2 dependence of the structure functions 
is included by an implementation of the DGLAP mechanism suitable for a cascade, so 
that the number of partons is not static, but varies in space and time as the collision 
of two nuclei evolves. The resulting parton distributions are presented and compared 
to distributions obtained with the VNI model. 



1 Talk given by V. Borchers in Quark Matter '99, May 13, 1999, Torino, Italy. For the transparencies, 
sec frittp: / / www.qm99.to. infn.it / program/qmprogram.html 



1 Introduction 



For a long time nuclear reactions have been described by phenomenological hadronic models 
quite successfully. At nuclear collision energies of the order of 10GeV-A it is already clear 
that partonic degrees of freedom will not be negligible; thus pure hadron cascades do not 
suffice in the SPS regime. 

Parton cascade models on the other hand provide a microscopic description of high energy 
collisions in terms of partons. The QCD-model of hard scatterings, the theoretical basis 
of parton cascades, limits them to reactions with high transferred momentum. The large 
number of involved particles (nucleons) in heavy ion collisions lets one assume that not all 
single nucleon-nucleon collisions will be hard. However, for sufficiently high energies of the 
incident nuclei one can expect that hard collisions are dominant. 

In the framework of a parton cascade non-perturbative effects can be either completely 
neglected or treated phenomenologically. For the model to be presented here we have chosen 
the first alternative, although this makes the model not applicable in the SPS energy regime. 

The first parton cascade model for heavy ion collisions, HIJING [fy], was formulated in 
momentum space only, and thus did not include a space-time description of the cascade 
evolution. A full phase-space description was first provided by Klaus Geiger's VNI [0, §], 
but this model is not manifestly Poincare- invariant. 

A description of the space-time evolution of the system is imperative if one aims to 
include effects due to the high density produced in the A — A collision. 

2 Parton cascades and Poincare- covariance 

The obvious way to enhance a parton cascade with a space-time description would be to 
choose an observer Lorentz frame and then demand that two particles i and j interact 
instantaneously as seen in this frame. It is clear that such a procedure makes the ordering 
of the interactions frame-dependent. This has been recognized as a problem of cascades for 
more than 15 years: it was found in hadronic cascades even at low energy (« 1 GeV-A) that 
the collision sequences and the resulting distributions depend noticeably on the reference 
frame (LAB or CMS) in which the simulation is run ||. 

The frame-dependence gets worse if one also makes the interaction criterion dependent 
on the spatial distance of particles []. 

If, on the other hand, we insist on Poincare-covariance we run into trouble because of 
the No Interaction Theorem ||, which asserts that interacting particles cannot be described 
Poincare-covariantly within a classical 6N-dimensional Hamiltonian theory of point particles. 

2.1 Poincare- Covariant Dynamics 

The solution of this problem chosen in the present work is the Poincare-Covariant Dynamics 
(PCD) approach described previously 0. In this approach the phase space is extended to 

2 It may be worthwhile to point out an additional frame dependence inherent in any description involving 
partons: the parton concept itself is frame-dependent, since it is well-dehned in the infinite- momentum frame 
only. 
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8N dimensions (N is the number of particles). Position and momentum vectors, r^(s) and 
Pi^(s) are parametrized by one Lorentz-invariant parameter s (which has no obvious physical 
meaning) . 

The particle motion and interactions are determined by the Lorentz-scalar 'Hamiltonian' 



N 



2 „2 



where V is a Lorentz-scalar pseudo-potential. The equations of motion are 

± ri (s) = {H, ri } = -?ji 

d . . r rr 1 dif 
_ Pi(s ) = {^ Pi } = + _ . 

Since the evolution proceeds in the invariant parameter s, every particle carries its own 
time ti = ri°(s). Note that in this formalism particles are classically off-shell when within 
the range of the quasi-potential V. 



3 Dynamics of the model 

Applying the concept of PCD to the Poincare Covariant Parton Cascade (PCPC), we simplify 
the dynamics by the following prescriptions: 

• Partons are free between collisions. 

• Binary interactions between partons i and j can occur when their 4-distance dij defined 

as 



dij ■= y - yr - —^rVJ w hh r := n - rj , p := p { + pj (1) 

has a minimum and will occur if this distance is smaller than the interaction distance 
^/er tot /7r. Note that d^ is an Lorentz-invariant way of writing the spatial (3-)distance 
in the center of momentum frame. 

Since the interaction criterion depends on Lorentz scalars only, the ordering of interactions 
is frame-independent. 



4 Initialization 

For setting up an initial state of the cascade, we define initial parton distributions in two 
steps, successively resolving the nuclei into nucleons and the nucleons into partons. First 
nucleonic Woods-Saxon distributions (in momentum and coordinate space) are used, and 



3 



the nucleons are boosted to the rapidities appropriate for the particular collision. Then the 
flavor and momentum of each parton are determined as follows: the parton momenta are 



p v = (y(xp/ ud ) 2 + P± 2 +/i 2 , P± , xp z 

The longitudinal momenta p z = xP z nncl of the partons are determined, together with their 
flavors, from the parton distribution functions f(x, Q 2 ) in the GRV94LO || parametrization. 
The distribution functions are evaluated at a fixed (and small) Qo , i.e. at a fixed scale. The 
parton transverse momenta \p±\ are Gauss-distributed (mean value 0, width 0.3 GeV). 

The one free quantity remaining is the parton's off-shell mass p, = \/m 2 — q 2 . This is 
fixed by the fact |/3 P | = |^ nucleon | _ This feature guarantees that partons originating 
from the same mother nucleon move together initially, thus modeling the constraint that 
initially partons are confined in the nucleons. Note how well PCD with its concept of 
potential-dependent masses is suited to describe such field-theoretic features in the classical 
framework of a cascade model. 

The parton spatial coordinates r« are distributed exponentially in the rest system of the 
nucleon. The boosted distributions are thus Lorentz-contracted automatically, and ad-hoc 
modifications of the spatial parton distributions (viz. "Distributed Lorentz Contraction" or 
"boost invariant sea parton distributions") that were proposed to enhance the interaction 
rate are not needed in this model. Due to the covariant formulation of our model they do 
not significantly influence the cascade evolution. 




5 Parton Interactions 

As stated in (Sec. |3|), partons i and j may interact at s — if their invariant distance dij(s) 
as defined in ([]]) has a minimum value. At this moment the interaction scale Q 2 (taken to be 
— (t+u)), and from that the relevant total cross section <Jij for the two partons is determined. 
Q 2 is not only the scale at which the running QCD coupling constant is evaluated but also 
decides if the scattering is "hard". Then the total cross section is determined from QCD in 
tree level approximation with inclusion of massive quarks. Only if then dy < v/Vy/V, the 
scattering actually occurs. 

Although our model is based essentially on binary hard scatterings, some non-perturba- 
tive effects are included in the case in which the interaction scale Qh of two partons is 
larger than the scale Qq 2 at which the partons were initially resolved. In that case the scale 
is adjusted via a "DGLAP scale evolution" |TI| [IT], |T^, [13], [T5|, This way 2 — > n events 



are effectively incorporated. In the picture below the scale evolution of parton a from the 
initial scale Qo 2 up to the interaction scale Qh 2 is depicted. The resolution scale is increased 
by successively radiated secondary partons, and it is the final parton b which is subjected to 
the binary scattering process^. 



3 Note the numbering: A larger scale means higher resolution, i.e. more partons. 
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6 Results 



-Q 2 



As numerical results we show (Fig. 2.) rapidity and transverse momentum distributions 
for central p-p and Au-Au collisions at RHIC energies (at \fs = 200GeVA). They were 
calculated with a Monte Carlo simulation based on the described model. Since the model 
in its present version does not include a hadronization mechanism, we do not compare our 
results with experimental data, but rather with a different parton cascade model. We have 



chosen the above-mentioned VNI code by Klaus Geiger [|14[], version 4.12 with hadronization 
disabled (partons only), and soft interactions switched off. 

In the overall features VNI and PCPC agree reasonably well. For p-p both models show 
a marked dip at midrapidity (as is to be expected). For Au-Au this dip is smeared out in our 
results, but, surprisingly, in VNI it remains essentially unchanged. Note also the significant 
differences in the transverse momentum distributions at low p±. 

In comparing these results, it should be noted that in VNI the hadron distributions 
obtained in runs with the hadronization mechanism switched on, are markedly different 
from the parton distributions shown here. 

7 Summary 

The Poincare-Covariant Parton Cascade Model (PCPC) as presented here is still incomplete. 
In this stage, it incorporates the following features: 

• new cascade-type microscopic simulation model of Nucleus- Nucleus Collisions at RHIC 
energies 

• basic quantities are partons (quarks and gluons) moving in 8N-dimensional phase space 
according to Poincare-Covariant Dynamics (PCD) 

• Interactions: hard parton-parton scatterings with cross sections determined from per- 
turbative QCD with a DGLAP scale evolution. 

What is left to do next is: 

• to include a hadronization mechanism 

• to add a thermodynamic analysis of results in space and time 

• to include shadowing effects, e.g. by using the EKS98 parametrization of parton dis- 



tribution functions [17 
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These are currently under development. 

The PCPC code [in C ++ ] is obtainable from 
http : //nix . physik . uni-bremen . de/ tkp. 
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Figure 1: Comparison of the PCPC and VNI parton distributions for central collisions at 
t/s = 200GeV-A : p-p (top panels) and Au-Au (bottom panels). The transverse momentum 



distributions 
the right. 
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are displayed on the left (cut: \y\ < 1), rapidity distributions 
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